Phase Separation in Binary Fluid 
Mixtures with Continuously Ramped 
Temperature 

By M. E. Cates 1 , J. Vollmer 2 , A. Wagner 1 ' 3 , D. Vollmer 1 ' 4 

1 School of Physics, University of Edinburgh, King 's Buildings, Edinburgh EH9 

3JZ, UK 

2 Max-Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, 

Germany 

3 Department of Physics, North Dakota State University, Fargo, ND 58105, USA 
4 Institute for Physical Chemistry, University of Mainz, Welder-Weg 11, 55099 

Mainz, Germany 

We consider the demixing of a binary fluid mixture, under gravity, which is steadily 
driven into a two phase region by slowly ramping the temperature. We assume, as a 
first approximation, that the system remains spatially isothermal, and examine the 
interplay of two competing nonlinearitics. One of these arises because the supersat- 
uration is greatest far from the meniscus, creating inversions of the density which 
can lead to fluid motion; although isothermal, this is somewhat like the Benard 
problem (a single-phase fluid heated from below). The other is the intrinsic dif- 
fusive instability which results either in nucleation or in spinodal decomposition 
at large supersaturations. Experimental results on a simple binary mixture show 
interesting oscillations in heat capacity and optical properties for a wide range of 
ramp parameters. We argue that these oscillations arise under conditions where 
both nonlinearities are important. 

Keywords: demixing, gravity, temperature ramp, binary fluids 

1. Introduction 

If a binary fluid, made of two species that are miscible at high temperature, is 
suddenly quenched into the two-phase region, it starts to demix. Depending on 
composition and quench depth, the mechanism for demixing and subsequent do- 
main growth is either spinodal decomposition (amplification of small compositional 
fluctuations, at first by diffusion and then by fluid motion driven by interfacial 
tension) or nucleation and diffusive growth of small droplets of one phase in the 
other (Bray 1994, 2000; Onuki 2002). Unless the two fluids have exactly matched 
densities, gravity eventually takes over, once the domain (or droplet) size becomes 
comparable to a suitably defined capillary length (Onuki 2002). The details of this 
gravitational stage are not fully understood and involve interesting new physics 
such as 'lane formation' (Chan & Goldberg 1975, Aarts et al. 2002). Nonetheless, 
it is observed that, once gravity does intervene, fluids separate relatively rapidly 
leading finally to a flat horizontal meniscus between phases. 
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In the natural world, and even in most laboratory settings, rapid temperature 
changes are a relatively rare occurrence. The opposite case of a very slow tempera- 
ture ramp, though equally idealised, is arguably closer to most everyday instances 
of phase separation, and it is certainly important to understand this limiting case. 
However, the physics is more complex, because the system is now subject to con- 
tinuous driving, as opposed to being displaced instantaneously from equilibrium, 
and allowed to relax back towards it. 

If the ramp rate dT/dt is small enough, one can hope that the time taken for 
heat to diffuse across the sample is small compared to other time scales of intcrestf. 
If so, this allows us to treat the system as isothermal at any instant; for simplicity we 
shall do this here, although in the experiments that motivate this work, the spatial 
gradients of temperature may not in fact be negligible. Neglecting all such gradi- 
ents reduces the problem from a double-diffusive one (Brandt & Fernando 1995), 
where diffusion of heat and composition arc both important, into one involving 
compositional diffusion only. 

Despite this the problem is still 'doubly nonlinear'. The first nonlinearity is 
standard, and arises from the coupling between density differences (caused by com- 
position deviations) to gravity. But the fact that the system can show phase sepa- 
ration requires that, even without this coupling, compositional diffusion is already 
a highly nonlinear process: a linearized diffusion equation obviously cannot yield 
spinodal decomposition and/or nuclcation of droplets. The interplay of these two 
nonlinearities, for the case of a slow isothermal temperature ramp in a binary fluid 
system undergoing phase separation, is addressed in this paper. 

2. Experimental motivation 

Vollmer et al. (2002) reported differential calorimctric studies of the phase separa- 
tion of a binary fluid mixture comprising water and 2 butoxy ethanol, also called 
C4E1, at low ramp rates (a few Kelvins per hour). The system has a lower conso- 
lute point and therefore demixes when heated up. The experiments revealed a quite 
unexpected effect: after an incubation period during which the phase boundary is 
first crossed and a meniscus forms, the observed heat capacity passes through a 
series of oscillations (typically about six in number) before eventually decaying to 
a smooth curve (Figure 1, upper left). The period and overall duration of the os- 
cillations depends on heating rate; it also depends on sample geometry, and on the 
composition of the initial state. 

An experiment in which heating was stopped for several hours midway through 
the oscillations (Figure 1, upper right) showed that these resumed when heating 
was recommenced, with only a small time delay; their period was not affected. 
Moreover, the eventual decay of the oscillations occurred exactly as though no 
interruption had occurred. These observations establish that there is little inertia 
involved in the process (since the interruption brings all fluid motion to a halt), 
and also suggest that the final decay of the oscillations is intrinsic, rather than 
stemming (for example) from a gradual loss of coherence in different parts of the 
system. (The recommencement of heating after a long break would presumably 
restore any required coherence.) 

f This time Th obeys A y ~ ht t h where A y is the smallest dimension of the sample cell (its 
thickness) and kt the thermal diffusivity. 
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Figure 1. Heat capacity measurements and videomicroscopy for mixtures of C4E1 and 
water, (top left) Apparent heat capacity of mixtures of various initial compositions [(a) 
0.18; (b) 0.25; (c) 0.35; (d) 0.39; (e) 0.43; (f) 0.44 % volume fractionC 4 Ei] that undergo a 
temperature ramp of 2 K/hour. (top right) Comparison of a scan (0.4 C4Eiby volume, 
0.8 K/hour ramp rate) with one where the heating was stopped after the first three 
oscillations and resumed at a later time (black), (bottom) Videomicroscopy sequence 
for a heating rate of 10 K/hour and an initial composition of 32 % C4E1 by volume. 
The temperatures (°C) are (a) 48.8, (b) 50.5, (c) 51.1, (d) 51.2, (e) 51.3, (f) 52.4, (g) 
53.4, (h) 57.9, (i) 58.8, (j) 59.5. The first oscillation in the sequence is visible in the final 
three frames. In frame (b) the sample is white. (The upper left and bottom figures are 
reproduced from Vollmer et al. (2002).) 
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By observing the upper phase present in these samples using contrast enhanced 
videomicroscopy (Figure 1, bottom) the oscillations were associated with the fol- 
lowing cycle. One observes an onset of high turbidity as large numbers of small 
droplets are suddenly created; these coarsen, and are soon cleared from the system 
by gravity. Bulk fluid motion is visible during this sedimentation stage; the system 
is then left in a quiet, transparent state with few droplets present. Soon after, the 
process repeats. In the quiet interludes between episodes of droplet formation, it is 
observed that fluid motion continues. 



3. Theoretical Considerations 

(a) Equations of motion 

We introduce a minimal model which is not strictly realizable in experiment 
(although it can more nearly be achieved in computer simulation; see Vollmer et ai, 
unpublished). This model addresses an isothermal, thermodynamically symmetric 
binary fluid, whose kinematic viscosity v (i.e., the diffusivity of momentum) is 
independent of composition and temperature, and whose mass density p depends 
on the composition variable (f> through a temperature independent parameter a — 
p~ 1 dp/d(j). The equilibrium phase diagram is symmetric about <ft — 4>c = 0, with the 
compositions of the two coexisting phases (at the binodal) obeying = ±(j> (T). 
The fluid is taken to be incompressible (V.u = 0). 

The Navier Stokes equation for the fluid velocity u may then be written 

^ + u.Vu + cm + u.V^j = z/V 2 u + agfoqiVz - V{p/p + gz) (3.1) 

Here g denotes gravity, acting towards — z; p is the fluid pressure; and if) — ((f> — 
<f>) /4>o is a relative supersaturation that measures the deviation of the local compo- 
sition from the equilibrium value (j> . This holds in regions of positive 0; otherwise 
the sign of <f> is reversed. (Note that a state of zero if> can describe two phases at 
equilibrium, with a sharp meniscus in between.) In what follows, we shall assume 
that the effects of composition on fluid density (proportional to a) mainly enter 
through the term in g on the right hand side: this term represents a buoyancy ef- 
fect in which local composition affects the body force acting directly on the fluid. 
The other place that a enters, on the left of the equation, represents the fact that 
changes in mass density caused by varying composition alter the acceleration of a 
fluid element subject to any particular set of forces. We shall neglect this below. 

The Navier Stokes equation is coupled to the nonlinear advection-diffusion equa- 
tion (c/. Appendix A for a derivation) 

^ + U .V0 = v(Df(4>)V4>) +(i-4>)£ (3.2) 

The comoving derivative on the left accounts for 'advection': the process whereby 
composition is transported by bulk fluid flow. The first term on the right represents 
diffusive currents, and the second is a source term for the supersaturation <f>. In 
contrast to the linear diffusion equations appearing in other settings, this equation 
has two remarkable features: 
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(i) The variable denotes 4>q 1 d(j)Q / dt and is a measure of the ramp rate. The 
source term, in which £ appears, characterizes the change of composition (j>o(T) 
with time during the temperature ramp, which gives rise to a constant increase 
of the relative supersaturation <j> under conditions where no diffusion or advection 
take place. As far as <j> is concerned, the effect of the temperature ramp is therefore 
to create supersaturation, at a rate controlled by £, throughout the sample. 

(ii) Equation 3.2 involves an effective nonlinear diffusivity Df(<j)), which decreases 
with the relative supersaturation <f), becoming negative for <j> > <f> 8 (Bray 1994). In 
particular D is chosen in such a way that /(0) = 1 and f((j> s ) — at the spinodal 
<p s where any state of uniform <f> > (f> s becomes locally unstable. Henceforth, we 
assume, for simplicity, that D and /(</>) do not depend directly on temperature, so 
that the only T-dependence of the diffusive dynamics enters parametrically through 
the definition, (p = 1 — (f)/4>o(T), of the relative supersaturation. 



(b) Dimensionless parameters 

The equations of motion can be made dimensionless by introducing suitable 
scales of mass, length and time. As a length scale we choose the height of the 
sample A = A z as the unit of length, and A 2 / D as time unit. A mass unit is pA 3 but 
interestingly, none is really needed. The fluid density in equation 3.1 factors through 
all terms in the equation, except for the one involving the pressure, V(p/p + gz), 
whose sole purpose is to maintain incompressibility. We will write this term as a 
dimensionless gradient Vp, but could equally well eliminate it by projecting the 
Navier Stokes equation onto incompressible flows in a standard fashion (see e.g. 
Onuki 2002). 

With these units, allowing for the various assumptions outlined above, and fur- 
thermore imposing a constant 'ramp rate' £, we find that the equations of motion 
take the following form: 



du. 

— + u.Vu = JVi 

at 



V 2 u + {Nq/N^Wz - Vp (3.3) 



^ + u.V0 - V (f($)V$) + (1 - $)N 2 (3.4) 



with three dimensionless parameter groups: 

N 1 = v/D ; N 2 =£A 2 /D ; ^ = agcj> h?/Dv (3.5) 

iV2 

Of these, N\ is a material parameter of the binary fluid (the ratio of momentum to 
particle diffusivities). N 2 can be thought of as a dimensionless ramp rate (although 
a different one, which does not depend on A, will be introduced later on). The group 
we have denoted N /N 2 is, in the same sense, a dimensionless gravity parameter. 
But we will find below that the product of this with N 2 plays a special role in the 
theory; that product deserves its own name, and we call it N = agcj) ^A 5 / D 2 v. 
Although these three groups can be combined in numerous ways to create new 
dimensionless numbers, they are sufficient to fully describe the parameter space in 
our idealized problem, for a given sample shape and given f ((/)). 

The statement that they are sufficient assumes that no important physics has 
been left out of the equations of motion 3.3 and 3.4. Perhaps the strongest candi- 
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dates for missing physics are (a) thermal gradients, which we have neglected from 
the outset; and (b) interfacial tension, which arise locally once the nonlinear regime 
of droplet formation is encountered, and are present at all times at the meniscus 
between phases. (Appendix A describes how a tension could be included.) 

One dangerous-looking combination of these neglected effects is when interfacial 
tension and thermal gradients combine to create Marangoni stresses (that is, spatial 
gradients of the tension); these are known to be implicated in several interfacial 
instabilities (Sternling & Scriven 1959, Davis 1987). While it would be imprudent 
to rule out an important role for Marangoni stresses and other thermal-gradient 
related effects, we do think it worth neglecting these in the first instance. 

(c) Relation to a Benard problem: Advective instability 

Consider first the limiting case where /(</>) = 1, or equivalently, where </> remains 
infinitesimal. Once the meniscus has formed, we can envisage a static diffusive state 
in which a steady current of supersaturation flows from the upper phase towards 
the lower and vice versa. (Note that, in such a state, <j> vanishes on the meniscus 
itself.) This current arises from the source term N2 in equation 3.4, which creates 
supersaturation uniformly through space. However, were such a current to come 
from a localized source that maintained constant supersaturation (f>* at the top of 
the sample, equations 3.3 and 3.4 would be isomorphic to those for the standard 
(Rayleigh-)Benard problem (see e.g., Faber 1995). The latter concerns a single phase 
fluid heated from below; this creates an inverted density gradient (with the denser 
fluid on top) so that above a certain heating rate the system becomes unstable. 
The role of the local temperature field, in the Benard problem, is played in ours 
by the supersaturation; that of the temperature difference between the plates, by 
<f>* . The correspondence applies whatever the sign of the expansion coefficient a: 
the less dense phase is always on top, so that supersaturation in this phase, which 
is largest far from the meniscus, always leads to a density inversion (and the same 
happens, symmetrically, in the lower phase too). 

The presence of a distributed source term, rather than driving at the bound- 
ary, gives an extra complication to this analogy due to the nontrivial form of the 
quiescent profile. This strongly alters the details of any stability analysis, but not 
the basic ideas. First, note that because the Benard problem is linear in the diffu- 
sive sector, it is governed by only two dimensionless groups, not three. This is the 
case in our problem also, so long as <fr (controlled by N2) remains small enough. 
The two parameters normally chosen for the Benard problem are the Prandtl num- 
ber, which is the direct analogue of our Aq, and the Rayleigh number, which in 
our language is Ra = {Nq/N2)4>*. The latter expression pertains to the 'standard' 
case of a constant supersaturation <p* at the upper boundary. When the source is 
distributed uniformly, this characteristic supersaturation scale instead depends on 
both the ramp rate and the system size, as 4>* ~ £A 2 /D = N2. This follows from the 
linearized diffusion equation: <j>* is the supersaturation that builds up at a distance 
of order A from the meniscus within the time it takes for diffusive relaxation over 
that distance. Combining these results, we may therefore identify No itself as the 
direct analogue of a Rayleigh number for our variant of the Benard problem. 

Much is known about the Benard problem, and using it we may now guess, with 
reasonable confidence, what will happen as the ramp rate £ is increased. The static 
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Figure 2. Suggested schematic phase diagram for an isothermally ramped binary fluid. 
The left panel gives the phase diagram in terms of the dimensionless control parameters 
A^o and N2, while the right is given in terms of A and £; both in log/log plot. The regions 
with qualitatively different behaviour are marked as S = stable, A = advectively unstable, 
D = diffusively unstable (perhaps with steady sedimentation), AD = advectively and 
diffusively unstable. Crossover lines from D to AD and from A to AD probably lie within 
the shaded regions shown. (In practice, of course, the boundary of the S region will not 
really have a sharp corner but will be a smooth curve, with the straight lines shown here 
as asymptotes. Likewise the 'crossover lines' need not meet the phase boundary precisely 
at this corner.) 

diffusive profile will be stable until Nq exceeds a critical value A^q (for the standard 
Benard problem, this is about 1700); thereafter, circulating fluid rolls will arise. 
This transition is represented in the phase diagram Figure 2 (left) as the transition 
from the region marked S to A. As A^ increases further, these rolls will undergo 
a series of bifurcations through various steady and perhaps unsteady (including 
oscillatory) states, with the details of this process controlled by N%. Eventually at 
some much larger No (also dependent on Ni), chaotic motion will set in, leading 
finally to fluid turbulence. 

Note in passing that, although the rolls in Benard are normally called 'convection 
rolls', they are properly 'advection rolls', in the following sense. The onset of these 
rolls is controlled by the advective nonlinearity (transport of supersaturation by 
fluid, which is the u.V0 term on the left side of equation 3.4) not the convective one 
(transport of momentum by fluid, which is the u.Vu in equation 3.3). For typical 
fluid parameters (Ni ~ 1 or more) one must move rather far along the sequence 
of transitions beyond N§ before the Reynolds number becomes large enough that 
convective, rather than advective, nonlinearity becomes important. 

(d) Onset of droplets: Diffusive instability 

All of the preceding subsection refers to the case where the characteristic super- 
saturation <p* remains small enough that the diffusion is linear. Let us now consider 
the opposite case where Nq remains small, but <j> does not. 
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For small N we expect no onset of advection rolls before the diffusive nonlin- 
earity kicks in. Within a linearized framework, we found <p* ~ N 2 , so the onset 
of the diffusive nonlincarity, if present on its own, must be governed simply by 
N 2 = £A 2 /D. When this reaches some critical value N 2 = N 2 (which in general 
depends on the shape of /(</>)), the stationary diffusion profile becomes unstable, 
which is indicated in the phase diagram (Figure 2, left) by the transition from the 
region S to D. In region D droplets are formed by either a spinodal or a nucleation 
process. The instability must occur at or before the point where the supersaturation 
in the static profile reaches (f> s ; this yields an upper bound N 2 < N^ (at which the 
region of highest supersaturation is locally unstable to spinodal decomposition). 
A particular choice of / yields N 2 — N^ 8 ~ 1 at small enough N (J Vollmer, 
unpublished), which appears to be adequate for our present purposes. In physical 
terms this result states that the diffusive instability kicks in whenever the height of 
the sample exceeds the length scale (D^) 1 / 2 beyond which diffusion can no longer 
compete with the homogeneous growth of supersaturation characterized by £. 

It is not completely clear what should happen once droplet formation begins. 
If gravity is strictly zero (N — 0) the system will coarsen without ever forming a 
horizontal meniscus; the presence of the ramp may cause domain formation within 
existing domains, reminiscent of patterns seen in certain binary systems undergo- 
ing simultaneous phase separation and polymerization (Clarke et al. 1995). (This 
process might conceivably lead to oscillatory states; the details certainly depend 
on N\.) In any case, once gravity is switched on, sedimentation of the nucleated 
domains will occur. By gathering supersaturation into localized droplets, nucleation 
is likely to result in the onset of advective nonlinearity sooner than would be the 
case without it; that is, before the line N = N§ is crossed (gray region in sector 
D). In the first instance the result could merely be steady sedimentation of the 
droplets; although this is, technically, an 'advective nonlinearity' (since it involves 
the u.V0 term in equation 3.4) it is relatively benign. By relaxing supersaturation, 
this could cause the threshold for bulk flow to move somewhat beyond N§; but 
such sedimentation, once present, is liable to give rise to fronts (Russel et al 1989) 
and other collective nonlinear behaviour (Chaikin 2000), distinct from those of the 
Benard problem. In summary, the onset of advective nonlinearity (beyond simple 
sedimentation), in a state where the diffusive one is already active, should occur at 
a threshold which is likely to lie somewhere between the diffusive instability line 
(N 2 = N 2 ) and the continuation of the advective instability line (No = N§). 

By the same token, starting with Nq > Nq and N 2 small, the presence of Benard 
rolls will clearly alter the criterion for onset of droplet formation. Here one can 
argue that the Benard rolls are more efficient than pure diffusion at transporting 
supersaturation, so that the onset of the diffusive instability will be delayed beyond 
the line N 2 = N 2 as was calculated for small No- On the other hand, for a given ramp 
rate £ (and fixed material parameters), we can find the cell size A* = A(Nq/Nq) 1 / 5 
for which the Benard instability would be just incipient. It would be surprising 
if, in a sample of size A > A*, the supersaturation far from the meniscus were 
actually less (for the same ramp rate and material parameters) than the value <fi* ~ 
7V2(A*/A) 2 pertaining to this smaller cell. This suggests that diffusive instability 
should set in at the latest upon achieving the condition N 2 (Nq/N ) 2 / 5 — N 2 , where 
the left side of the equation is the value of N 2 in a fictitious sample of height 
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A*. In summary, the onset of diffusive nonlincarity, in a state where the advective 
one is already active, should occur somewhere in the shaded region marking the 
transition from the A to the AD regions in the left diagram of Figure 2, whose 
boundaries are the continuation of the diffusive instability line (N 2 = N 2 ) and the 
line N 2 /N% ~ (N Q /N§) 2 / 5 . 



(e) Nonequilibrium phase diagram for experimental control parameters 

The above arguments, which were summarized in the left panel of Figure 2, 
predict a phase boundary between a stable diffusive regime (S) and an advectively 
unstable one (A) at No = N§; and a boundary between S and a diffusively unstable 
regime (D) at N 2 — N 2 . For comparison with experiment it is helpful to represent 
these lines on a plot that involves two more accessible parameters, the system size 
A and ramp rate £, as the horizontal and vertical axes (Figure 2, right). To this 
end we observe that the combinations N /N 2 = (A/A ) 3 and N^/Nq = (£/£o) 3 
only depend on material constants. In these expressions we have identified 'natural' 
units of length and heating rate A 3 , = Dv/ag4> and = a 2 g 2 '4>lD / V 2 = (Z5/A 2 ,) 3 , 
respectively, only depend on material constants. Consequently, on logarithmic scales 
the phase diagrams can be related by the linear transformation 



v-i §A'»s«u (3 ' 6) 

where A c and £ c are chosen so that the two boundaries of the stable region S intersect 
at (1, 1), i.e., the sharp corner visible in the right panel of Figure 2. Consequently, 
the S/A and S/D lines now meet at the point (A c , £ c ), where A c /A = (A^/A^f) -1 / 3 , 
and £ c /£ = (7V 2 C ) 5 / 3 (7V C )- 2 / 3 as defined previously. 

The estimates for the boundaries between the stable domain S and the advec- 
tively (A) and diffusively (D) unstable ones translate in this phase diagram into 
straight lines with slope -5 and -2, respectively. The limiting estimates made above 
for where regimes A and D should each cross into AD (the regime in which ad- 
vective and diffusive nonlinearities arc simultaneously strong) are shown by dashed 
lines as previously, and the uncertainties for these regions are again indicated by 
gray areas. The estimate of the upper bound for the appearance of nucleation in the 

2/5 

advective regime N 2 ~ N ' translates into a condition on the ramp rate £, and is 
independent of the sample height A (although it will, like other critical parameters, 
in general depend on the shape of the sample cell, and on /(</>)). 

Not shown in either phase diagram are additional, increasingly vague crossovers 
from A, D and AD into further regimes where the convective nonlinearity also 
becomes important. 



4. Oscillation Mechanism 

In the experiments showing oscillatory demixing in C4E1, that were outlined in §2, 
one observes fluid flow (without obvious signs of turbulence) and droplet formation 
simultaneously. This means that advective and diffusive nonlinearities are both 
involved, making it plausible that these experiments lie within the AD region of 
the phase diagram. Preliminary parameter estimates (Vollmer et a/., unpublished) 
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suggest that they may lie in a part of the AD region well to the right of, but not 
much above, the point (A c ,£ c ) on Figure 2 (right). However, in the experiment it 
is dT/dt that is held constant, not £, which can therefore drift slowly during the 
ramp, as can a, u, etc. These drifts may cause the experimental parameters to move 
gradually towards more stable values, causing the oscillations finally to cease. 

These remarks do not explain why the oscillations are there in the first place. 
Elsewhere we propose that in this part of the AD region, advection rolls form soon 
after creation of the meniscus, while the supersaturation <fi is still small. These rolls 
do not advect enough flux to maintain <fi below the threshold for droplet nucleation. 
Hence droplets appear; they then grow for a time, but before long, sedimentation 
sets in. The falling droplets soak up the remaining supersaturation, restoring the 
system to a state of small 4>, so that the cycle is ready to begin again. In Vollmer 
et al. (unpublished) we find this scenario consistent with preliminary computer 
simulation data, and broadly in line with experimental data on various systems. 

This explanation has something in common with a proposed mechanism for 
rhythmic deposition of precipitate in a convecting magma chamber cooled from 
above (Sparks et al. 1993), a problem that seems, however, to involve a significantly 
larger parameter space than the model considered here. 

5. Conclusion 

Even within the idealized model presented here, the interplay of gravity-driven 
advection and nonlinear diffusion presents an interesting theoretical challenge for 
the description of binary fluid demixing in the presence of a temperature ramp. 
The experimentally observed oscillations (Vollmer et al. 2002) appear in a regime 
where both nonlinearities are simultaneously large; this makes analysis difficult. 
Detailed computer simulations, now planned, may be a better way of finding out 
whether the physical ingredients retained in the idealized model (equations 3.3 and 
3.4) arc sufficient to explain the observed oscillations, or whether factors that we 
have omitted, such as Marangoni stresses arising from thermal gradients, play an 
essential role in this problem. 

J.V. acknowledges financial support from the Schloessmann Foundation of the Max Planck 
Society, and thanks Howard Stone and Michael Brenner for useful discussions. D. Vollmer 
thanks the EC for a Marie Curie Fellowship. 

Appendix A. Advect ion-diffusion equation for the 

supersaturation <p 

The evolution of the composition <f> of the mixture is governed by the advection- 
diffusion equation 

^ + u.V0 = -V J (Al) 

where the comoving derivative on the left accounts for the advection of composition 
by the bulk fluid flow; and the divergence of the diffusive current J<f, on the right 
hand side accounts for the diffusive decay of composition gradients. The current 
Jcf, = — MV/i, with M a mobility, is driven by the gradient of chemical potential n, 
which — in local equilibrium — obeys [i = 5F/5(j), where a suitable choice for the 
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free energy functional F[<j>] is 

F[cf>} = J dV [F o (0) + «(V^) 2 ] (A 2) 

This involves the free energy density F ((j)) describing homogeneous phases in equi- 
librium, and an energetic penalty k(V(/>) 2 for steep changes of composition. 

For nonvanishing k, the free energy functional A 2 supports static interfaces 
between phases which have small but finite width, and gives an equilibrium inter- 
facial tension that depends on k and F n ((ft>) (Bray 1994, Onuki 2002). Although 
intcrfacial tension can thus be included without much modification into the equa- 
tions of motion, this would introduce a fourth dimensionless control parameter into 
our model. We assume here that it plays no major role and omit the corresponding 
terms from the equations of motion, although in our preferred simulation algorithm 
for this type of problem (lattice Boltzmann; see Kendon 2001), they are actually 
incorporated in precisely the above way. 

Evaluating Al with n — and = —MV(SF/5(f>), introducing the abbre- 
viations D = M [^ 2 ^oM/<50 2 ] 0=0o , f{4>) = (M/D) 8 2 F [cf>]/5<i) 2 , and taking into 
account the temperature dependence of the factor <fi appearing in the dimensionless 
supersaturation <f> one obtains equation 3.2 for the evolution of qb. 

References 

Aarts, D. G. A. L.; van der Wiel, J. H. & Lekkerkerker, H. N. W. 2002 Interfacial dynamics 
and the static profile near a single wall in a model polymer mixture. J. Phys. Cond. 
Matt ** *** *** 

Brandt, A. & Fernando, H. J. S. (eds) 1995 Double diffusive convection. AGU Geophysics 

Monograph Series, Washington DC. 
Bray, A. J. 1994 Theory of phase-ordering kinetics. Adv. Phys. 43 357-459. 
Bray, A. J. 2000 Coarsening dynamics of nonequilibrium phase transitions. In Soft and 

fragile matter, nonequilibrium dynamics, metastability and flow (eds Cates, M. E. & 

Evans, M. R.) IOP Publishing, Bristol. 
Chaikin, P. M. 2000 Thermodynamics and hydrodynamics of hard spheres: the role of 

gravity. In Soft and fragile matter, nonequilibrium, dynamics, metastability and flow 

(eds Cates, M. E. & Evans, M. R.) IOP Publishing, Bristol. 
Chan, C. K. & Goldberg, W. I. 1987 Late-stage phase separation and fluid flow in a binary 

liquid mixture. Phys. Rev. Lett 58, 674-677. 
Clarke, N., McLeish, T. C. B & Jenkins, S. D. 1995 Phase-behavior of linear branched 

polymer blends. Macromolecules 28, 4650-4659. 
Davis, S. H. 1987 Thermocapillary instabilities. Ann. Rev. Fluid Mech. 19, 403-435. 
Faber, T. E. 1995 Fluid mechanics for physicists, Cambridge University Press, Cambridge. 
Kendon, V. M.; Cates, M. E., Pagonabarraga, I., Desplat, J.-C. & Bladon, P. 2001 Inertial 

effects in three-dimensional spinodal decomposition of a symmetric binary fluid mixture: 

a lattice Boltzmann study. J. Fluid Mech. 440, 147-203. 
Onuki, A. 2002 Phase transition dynamics, Cambridge University Press, Cambridge. 
Russel, W.B., Saville, DA. & Schowalter, W.R. 1989 Colloidal Dispersions, Chapter 12, 

Cambridge University Press, Cambridge. 
Sternling, C. V. & Scriven, L. E. 1959 Interfacial turbulence: Hydrodynamic instability 

and the Marangoni effect. A.I.Ch.E. J. 5, 514-523. 



Article submitted to Royal Society 



12 



M. E. Cates, J. Vollmer, A. Wagner, D. Vollmer 



Sparks, R. S., Huppert, H. E., Koyaguchi, T. & Hallworth, M. A. 1993 Origin of modal and 
rhythmic igneous layering by sedimentation in a convecting magma chamber. Nature 
61, 246-249. 

Vollmer, D., Vollmer J. & Wagner A. J. 2002 Oscillatory kinetics of phase separation in a 
binary mixture under constant heating. Phys. Chem. Chem. Phys. 4, 1380-1385. 



Article submitted to Royal Society 



